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We present an analysis of the model by Gnanaa 
[1999] for the pycnocline depth in the ocean. An analytic 
solution for the overturning strength as a function of the 
meridional pressure difference is derived and used to discuss 
their mutual scaling. We show that scaling occurs only in 
two unphysical regimes of the model. In the absence of the 
Southern Ocean (SO) processes, i.e. for a northern overturn- 
ing cell, the volume transport is proportional to the square 
root of the pressure difference. Linear scaling is seen when 
the overturning is restricted entirely to the SO, i.e. when 
no northern downwelling exists. For comparison, we present 
simulations with the coupled climate model CLIMBER-3a 
which show linear scaling over a large regime of pressure dif- 
ferences in the North Atlantic (NA). We conclude that the 
pycnocline model is not able to reproduce the linear scal- 
ing between its two central variables, pressure and volume 
transport. 



1. Introduction 

The meridional overturning circulation in the Atlantic is 
a central challenge to our understanding of global climate 
dynamics. Gnanadesikan [1999] {G99 hereafter) presented 
a model for the deep meridional circulation in terms of the 
pycnocline depth (PD). This idealised model has been under 
intense investigation as a possible paradigm for the merid- 
ional overturning circulation [Gnanadesikan and Hallberg, 
2000; Saenko and Weaver, 2002; Gnanadesikan et al., 2002]. 
Furthermore it has been used to investigate the qualitative 
importance of different physical feedbacks on the oceanic 
circulation [Klmger et al., 2003, Gnanadesikan et al., 2003; 
Kamenkovich and Sarachik, 2004]. A qualitative feature 
of the deep meridional overturning circulation is the scal- 
ing relation between the volume transport and the merid- 
ional density difference in the Atlantic [Bryan 1987]. Pick- 
ing up Bryan's scaling arguments but assuming a constant 
PD in the Atlantic Rahmstorf [1996] proposed a linear re- 
lation which he demonstrated in the oceanic general cir- 
culation model (GCM) MOM-2. Park [1999] and Scott et 
al. [1999] derived the same scaling in a Stommel-type box 
model. GCM simulations of the ocean suggest that this lin- 
ear relation carries over from the density to the pressure 
difference [Hughes and Weaver, 1994; Thorpe et al., 2001]. 
The linear scaling relation between pressure difference and 
maximum overturning strength has since been demonstrated 
to be a robust feature in oceanic GCM simulations ( Griesel, 
pers. comm., 2004). In sec. 4 we present simulations with 
the coupled climate model CLIMBER-3a further supporting 
these findings. 
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The G99 model contains four physical processes which 
influence the PD in the ocean. The balance of the pres- 
sure gradient in the North Atlantic and the frictional forces 
within the boundary currents leads to an equation for the 
northward volume transport 
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The pressure gradient is parameterised through the density 
difference in the NA Ap, the north-south distance Ly 
which the gradient occurs and the PD D. 



Ap = gDAp 



(2) 



The constant 7„ combines lJ,"' with f}, p and C (the merid- 
ional derivative of the Coriolis parameter /, the density and 
a proportionality constant of order one), g is the gravity 
constant. The quadratic dependence on D occurs due to 
the vertical integration in order to obtain a volume trans- 
port. In the SO the model includes the Drake passage effect 
through a wind-driven upwelling which does not explicitly 
depend on the PD ri*^' — {Lxt) / {pf) = 27e. r and Lj, 
are the wind stress in the SO and the circumference around 
Earth at the latitude of Drake Passage. Additionally G99 
includes an eddy induced return flow 
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where v^d is the transport velocity which G99 parameterised 
following Gent and McWtlUams [1990] while we focus here 
on its dependence on the PD. The fourth term in the model 
is associated with low-latitudinal upwelling described by an 
advection-diffusion balance wdzP = K^dzzP in the tropics 
which yields 
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where and Au are the diapycnal diffusivity and the hori- 
zontal area of upwelling, respectively. All non-negative con- 
stants 7^ have been introduced for convenience. Note that 
the underlying assumption of the model is that these four 
process can be described using the same value D for the 
PD throughout the Atlantic. Equ. (1) requires futhermore 
that the vertical extension of the northward volume flow is 
also given by D. Accepting these assumptions, the conser- 
vation of volume then results in the governing equation of 
the model 



= jngAp ■ D + ■jgm ■ D - 27e ■ -D - 7u 



(5) 



It can be shown that for all parameter settings the model 
has at most one solution with non-negative PD. In sec. 2 
we give this solution analytically in terms of the volume 
transport T„ as a function of the pressure difference Ap and 
discuss, in sec. 3, its scaling with Ap. In sec. 4 we compare 
the results with simulations with the coupled climate model 
CLIMBER-3a. 
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2. Solution for the Volume Transport T„ 

In order to obtain an analytic solution of the model we 
rewrite equ. (1) to get an expression of the PD as a function 
of volume transport T„ and pressure difference Ap 

D = T„/(7„Ap) (6) 

In the most interesting case of non-zero volume transport, 
T„ ^ 0, we can insert the equality (6) into the volume con- 
servation equ. (5) to get 

= r„-F^— -27e-7^7„7^ 7 

7n Ap r„ 

Multiplying by r„ yields a quadratic equation in T„ with 
two solutions of which only one is non-negative 

^" = ^ ll^I^ (> + + In (rtn^P + 7.m)) (8) 

Note that despite the fact that the governing equ. (5) is cu- 
bic in D, the model does have at most one physical solution 
given by equ. (8). The model does therefore not bear the 
possibility of multiple stable modes of the deep meridional 
overturning circulation as suggested by simulations with cli- 
mate models of different complexity [Stommel, 1961; Man- 
abe and Stouffer, 1988; Rahrnstorf, 1995, 1996; Ganopolski 
et.al., 2001; Prange et al., 2003]. This is to be expected given 
that the model does not include a salt-advection feedback as 
proposed by Stommel [1961]. Fig. 1 shows the solution for 
different diapycnal mixing coefficients K,, . The results were 
obtained using the numerical values given by G99. Note 
that the solution (8) depends continuously on the diapycnal 
mixing coefficient oc 7^. No change in the quality of 
the solution (8) occurs in the absence of the low-latitudinal 



upwelling, where 

j,(K„=0) ^ 27e7nAp ^ 7nAp _ ^(e) 
7nAp-|-7gm 'yn^P + '^gm 

In contrast to the behaviour for vanishing K^, the elimina- 
tion of the SO processes changes the quality of the solution 
as can be seen from equ. (8) and will be discussed in the 
next section. 

3. Scaling of the Volume Transport T„ 

Next, let us discuss the scaling of the volume transport 

Tn with the meridional pressure difference Ap. First, con- 
sider the situation without the SO processes, i.e. ri*^^ = 
rp(am) _ Q rpj^g scaling can be obtained from the general 
solution in equ. (8) with 7^ = 7gm = 0. More illustrative 
is the derivation from the original equations for the volume 
transport (1) and (4). The fact that the northern down- 
welling has to be balanced by the low-latitude upwelling 
Tn = 'ynD ■ Ap = Tu = "iu/D implies that 

i.e. the PD decreases with increasing pressure difference in 
the NA. Using this expression to replace D in the parame- 
terisation of T„ in equ. (1) yields 

ri""^") = V7.7n • Ap ~ ^ (11) 

In connection with equ. (2) we get the scaling ~ 

(Ap)^/^ which was derived first by Bryan [1987]. 

Next let us add the SO winds, but neglect the eddy-induced 

return flow, i.e. 7gm = 0. The solution (8) then be- 
comes Tn = + + 7u7nAp, whlch goes to a constant 
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Figure 1. The analytic solution of the conceptual model 
for the volume transport Tn as a function of the merid- 
ional pressure difference Ap for different values of the 
diapycnal diffusivity Kv The dots represent simulations 
with the coupled model CLIMBER-3a which contains an 
oceanic GCM. The simulations exhibit a linear scaling in 
contrast to the conceptual model. 
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Figure 2. The same curves as in fig. 1 in double log- 
arithmic scale. The dashed line corresponds to the so- 
lution without the SO processes and shows a one-half 
scaling Tn ~ VAp. In contrast to this the solutions 
which include the SO processes show a linear relation- 
ship Tn ~ Ap. 
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T-a — Ts'^ = const, in the wind-driven limit, i.e. for 

small vertical difltusivity 7^ <C 7e/(7nAp). As expected no 
scaling between pressure difference and volume transport is 
observed in this case. 

In order to discuss the scaling behaviour in the presence 
of both SO processes in the model, we plot the solutions 
shown in fig. 1 in double logarithmic scale in fig. 2. For 
small Ap all solutions which include the SO processes have 
slope one which corresponds to a linear scaling of the vol- 
ume transport with the pressure difference. For comparison 
the solution without the SO processes from equ. (11) has 
been included as the solid curve in fig. 2 showing the one- 
half slope. This result can be understood from the general 
solution in equation (8) which also sets the scale Aps for 
which the linear relation holds. For Ap <§; Aps = 'Ygm/^n 
the solution (8) becomes 



T„ = 



1 + 



7e 



• Ap ~ Ap, (12) 



i.e. Tn is linear in the meridional pressure difference Ap. 
Using the numerical values given by G99, we obtain an es- 
timate for the pressure scale Aps = 31.25 hPa which is 
consistent with the scaling seen in fig. 2. Simulations with 

the oceanic general circulation model MOM-3 show a linear 
scaling of T„ with Ap for a variety of parameter settings, 
including the case of zero diapycnal mixing (Griesel, (pcrs. 
comm., 2004). The pressure scale in these simulations is of 
the order of Aps « 50hPa which is in good agreement with 
the above estimate. 

The physical meaning of the scaling regime is seen when 
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multiplying Aps with D which gives T„ <C 
means that the scaling occurs only when the circulation is 

completely dominated by the SO processes, i.e. when the 
eddy-induced return flow in the Southern Ocean is nmch 
stronger than the downwelling in the NA. This situation 
is not consistent with the underlying physical assumption of 
the model of an interhemispheric meridional overturning cir- 
culation and it does not describe the observed circulation in 
the ocean. From equ. (6) and (12) we can see that in the lin- 
ear scaling regime the pycnocline depth does not vary with 
the pressure difference, in contrast to the situation without 
SO processes (equ. (10)) where D decreases with Ap. From 
equ. (2) we see that for constant D the pressure difference 
scales in the same way as the density difference Ap ^ Ap 
making the linear scaling a simple consequence of the initial 
assumption that T„ oc Ap (equ. 1). 



4. Comparison with simulations 

The linear relationship between the maximum overturn- 
ing strength and the density difference Ap which was ob- 
served by Rahmstorf [1996] in a oceanic GCM is reflected 
in the parameterization of the northern downwelling in 
equ. (1). In the conceptual model, however, it does not 
carry over to the pressure difference, as was shown in the 
previous section. In order to check this scaling we carried 
out sinmlations with the coupled climate model CLIMBER- 
3a. The model contains an atmosphere and a sea-ice mod- 
ule as well as the oceanic general circulation model MOM- 
3. The effect of baroclinic eddies was included through 
a parameterization following Gent and McWilliams [1990] 
with a coefficient of Kgm = 2.5 • 10^cm^s~^. For a full de- 
scription of the model see Montoya et al. [2004]. Starting 
from the present day equilibrium simulation with a max- 
imum overturning strength of 12 Sv, we apply a negative 
salinity forcing of different strength to the NA convection 



sites (between 50°N and 80°N) as described in Levermann 
et al. [2004] . This leads to a decrease in the meridional pres- 
sure difference in the NA and therefore a weakening of the 
meridional overturning. A positive salinity forcing strength- 
ens the overturning and increases the pressure difference. 
Fig. 1 shows the simulations as black dots. The pressure 
was taken at a depth of 1500 m corresponding to the center 
of the overturning cell in the simulations. The differences 
were taken between the zonal average between 50° N and 
80°N and the zonal average between 20°N and 30°N. This 
corresponds with the meridional pressure difference in the 
NA that enters equ. (1). As seen in fig. 1 the maximum 
meridional overturning in the Atlantic scales linearly with 
the pressure difference in the NA in the simulations. The 
vertical diffusivity in coupled model was kept constant at 
Kv = 0.1 cm^s~^. Thus the simulations correspond to the 
dashed solution curve in fig. 1. Simulations and concep- 
tual model do neither agree quantitatively using the values 
suggested by G99 nor is the qualitative behaviour of the 
two main quantities (pressure and volume transport) repro- 
duced in the conceptual model. These results are supported 
by recent findings by Griesel (pers. comm., 2004) with an 
oceanic GCM. Their work shows that the linear scaling be- 
tween pressure and overturning strength is a robust feature. 
It is independent of changes to various parameters including 
the Gent and McWilliams diffusivity coefficients. In order 
to emphasize that fact that the linear scaling r„ ~ Ap cor- 
responds to constant D we plot in fig. 3 the PD as defined 
in G99 for our simulations. In contrast to G99 in an OGCM 
we find in our coupled model no significant variation of the 
PD for varying pressure difference. 



5. Conclusions 

By giving an analytic expression for the meridional over- 
turning strength Tn as a function of the meridional pres- 
sure difference Ap, we discuss the scaling of the two main 
quantities of the conceptual model introduced by G99. The 
model exhibits two scaling regimes which both correspond 
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Figure 3. Pycnocline depth for the simulations shown in 
fig. 1 as a function of the meridional pressure difference 
Ap in the NA. Definition and displayed depth range were 
taken as in G99. 
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to unphysical situations. Linear scaling occurs in a situa^ 
tion where the eddy-induced return flow is much stronger 
than the northern downwelling. This corresponds to a cir- 
culation which is localized entirely in the SO and in which 
all downward volume transport is due to the eddy-induced 
return flow. This situation is inconsistent with the physical 
assumption of an interliemispheric overturning cell underly- 
ing the model and the isopycnal nature of the return flow. 

The second scaling regime corresponds to a purely north- 
ern cell where the upwelling takes place entirely in low 
latitudes, described by an advection-diffusion balance. In 
this case the overturning is proportional to the square root 
of the pressure as reported by Bryan [1987]. The scaling 
was checked using the coupled climate model CLIMBER- 
3a with a parameter setup comparable to the conceptual 
model, i.e. including effects of baroclinic eddies following 
Gent and McWilliams [1990] and a vertical diffusivity of 
Kh = O.lcm^s"^. The simulations exhibit a linear scaling 
and therefore support previous studies [Hughes and Weaver, 
1994; Rahmstorf, 1996; Thorpe et al, 2001] with comprehen- 
sive climate models. The PD does not vary significantly as 
a function of the pressure difference in our simulations. 

We conclude that the conceptual model of the PD can 
not reproduce the scaling between its central variables, the 
pressure and the volume transport. Besides possible crit- 
icism regarding the specific parameterizations of the four 
physical processes contained in the model, the assumption 
of a universal D for all these processes seems questionable. 
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